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Heterogeneous nucleation of a new bulk phase on a flat substrate can be associated with the 
surface phase transition called wetting transition. When this bulk heterogeneous nucleation occurs 
on a completely-wettable flat substrate with a zero contact angle, the classical nucleation theory 
predicts that the free energy barrier of nucleation vanishes. In fact, there always exist a critical 
nucleus and a free energy barrier as the first-order pre-wetting transition will occur even when the 
contact angle is zero. Furthermore, the critical nucleus changes its character from the critical nucleus 
of surface phase transition below bulk coexistence (undersaturation) to the critical nucleus of bulk 
heterogeneous nucleation above the coexistence (oversaturation) when it crosses the coexistence. 
Recently, Sear [J.Chem.Phys 129, 164510 (2008)] has shown by a direct numerical calculation of 
nucleation rate that the nucleus does not notice this change when it crosses the coexistence. In 
our work the morphology and the work of formation of critical nucleus on a completely-wettable 
substrate are re-examined across the coexistence using the interface-displacement model. Indeed, 
the morphology and the work of formation changes continuously at the coexistence. Our results 
support the prediction of Sear and will rekindle the interest on heterogeneous nucleation on a 
completely-wettable substrate. 



PACS numbers: 64.60.-i, 64.60.1 
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I. INTRODUCTION 

The nucleation which occurs within the bulk 
metastable material is called homogeneous nucleation^. 
However, it is frequently assisted by the presence of a sur- 
face. The nucleation in this case is called heterogeneous 



where 



deatior 



In the course of vapor to liquid nucleation, 



for example, a liquid droplet of semi-spherical shape will 
be formed on the substrate which attracts liquid within 
the oversaturated vapor. The shape of this critical nu- 
cleus of the liquid is characterized by the apparent con- 
tact angle 9 a shown in Fig. 1. 

The nucleation is characterized by the nucleation rate 
J, which is the number of critical nuclei formed per unit 
time per unit volume. Usually it is written in Arrhenius 
form 



J = A cxp 



W 



(1) 



where A is a kinetic pre-exponential factor which is be- 
lieved to be weakly dependent on temperature T, fee is 
the Boltzmann's constant, and W is the reversible work 
of formation of critical nucleus 

According to the classical nucleation theory (CNT)£, 
the work of formation (free energy barrier) of the hetero- 
geneous nucleation Wietero is expressed using the appar- 
ent contact angle 9 a and the free energy barrier Wiomo 
of the homogeneous nucleation as 



/ (9 a ) = {9a - cos 9 a sin 9 a ) /it 



(3) 



for the two-dimensional semi-cylindrical nucleus, and 

f(9 a ) = (l-cos0 Q ) 2 (2 + costf a )/4. (4) 

for the three-dimensional axi-symmetric semi-spherical 
nucleus. The contact angle 9 a characterizes the inter- 
action between the liquid and the substrate. When 
7T > 9 a > the liquid is said to incompletely wet the 
substrate. When 9 a — 0, the liquid is said to completely 
wet the substrate. Since < / (9 a ) < 1, the presence of 
the substrate enhances the nucleation rate (Eq. ([1])) when 
the liquid incompletely wet the substrate (tt > 9 a > 0) 
as Whetero < Whomo- When the liquid completely wet 
the substrate (9 a = 0), CNT predicts that no nucleation 
barrier exist (Whetero = 0) as / (9 a = 0) = 0. 



Droplet 




Whetero = Whomo/ (9 a ) , 



(2) 



Substrate 



FIG. 1: An ideal critical nucleus on a substrate. The appar- 
ent contact angle 6 a is defined as the angle of the intersection 
of the cylindrical or spherical droplet surface with the effec- 
tive radius R e fi and the height L from the substrata and the 
wetting layer of the thickness l e . 
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This heterogeneous nucleation of a new bulk phase on 
a flat substrate can be associated with the surface phase 
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transition called wetting transition as both phase tran- 
sitions can occur on the same surface^ -. In fact, there 
always exists a critical nucleus and a free energy bar- 
rier of the first-order so-called prewetting transitioniir— 
even on a completely-wettable substrate. This transition 
occurs in the undersaturated vapor and separates states 
of thin and thick liquid films on a completely-wettable 
substrate with the zero contact angle. The shape of the 
critical nucleus is flat and pancakes-like^, and cannot be 
characterized simply by a contact angle. Therefore, in 
contrast to the prediction of CNT, the critical nucleus 
and the free energy barrier for the heterogeneous nucle- 
ation does exist even when the substrate is completely- 
wettable and the contact angle 6 a = (Eq. ([2])). By 
varying the chemical potential or the vapor pressure from 
undcrsaturation to oversaturation, the droplet changes 
its character from the critical nucleus of the prewetting 
surface phase transition to that of the vapor to liquid 
bulk phase transition. Recently, Sear— has shown by a 
direct numerical calculation of nucleation rate that the 
nucleus " does not notice" this change when it crosses the 
coexistence. 



In our work, the morphology and the work of formation 
of a critical nucleus rather than the nucleation ratei^ on 
a completely-wettable substrate are re-examined across 
the coexistence using the interface-displacement model 
(IDM) which has been successfully used for studying wet- 
ting phenomenal— , line tension o 16 ' 17 and even layering 
transitions^. In fact, various properties of critical nu- 
cleus of wetting layer rather than the bulk critical nu- 
cleus have already been studied. For example, the sta- 
bility and the critical exponents of the critical nucleus of 
wetting layer near the wetting and prewetting transitions 
have been studied using IDMi2r— and Nakanishi-Fishcr 
mode l 23 ' 24 . However those authors payed most attention 
to the nucleation of wetting transition in the undersat- 
urated vapor below the bulk coexistence. In fact, the 
heterogeneous nucleation of bulk phase transition occurs 
in the oversaturated vapor above the bulk coexistence. 
Therefore, we will pay attention to those properties which 
is relevant to the bulk heterogeneous nucleation, which 
only a few number of authors such as Talanqucr and Ox- 
tobyiS and Sear— have considered 



In Sec. II, we shall briefly recall the description the 
first-order surface phase transitions within the framework 
of IDM. In Sec. Ill we shall study the morphology and 
the thermodynamics of critical droplets. In particular, 
we shall focus our attention not to the scaling properties 
but on the relation between the heterogeneous nucleation 
and the prewetting transition. Section IV is devoted to 
conclusion. 



II. A SHORT REVIEW OF THE WETTING 
TRANSITION AND THE 
INTERFACE-DISPLACEMENT MODEL 

Within the interface displacement model (IDM) in d- 
dimensional apace the free energy of a fluid film of local 
thickness I (x) is given by^ ~ 21 ' 25 ~ — 



[I] = J 



7 



1 + (V/)' 



1/2 



-1 + V(Z)-/xZ 



(5) 

where 7 is the liquid-vapor surface tension, V (I) is the ef- 
fective interface potentia l 11 ' 13 ' 30 from the substrate, and 
(i denotes the deviation of the chemical potential from 
liquid-vapor coexistence such that fx = is the bulk co- 
existence and for [i > the vapor is oversaturated. In 
Eq. ((5j) we keep the nonlinear dependence on VZ. 

Sometimes it is useful to consider the full potential (j>(l) 
defined by 



0(Z) = V(Z)-/iZ. 



(0) 



instead of the effective interface potential V (I) and the 
chemical potential contribution fxl separately. This full 
potential </>(/) depends not only on the chemical potential 
\x but also on the temperature T. Figure [5] shows typical 
shapes of the full potential 4>{l) of a completely-wettable 
substrate. The full potential (f>(l) exhibits double-well 
shape typical to the first-order surface phase transition, 
and its two minima at l e and at L e correspond to the 
metastable thin and the stable thick wetting films. Fig- 
ure [2] indicates that the thick film become infinitely thick 
(L e = 00) when \x > 0. 

The schematic surface phase diagram is shown in Fig. [3] 
in the T-fi plan e 12 ' 13 ' 16 ' 17 ' 19 ' 21 . Since the local minimum 
V (Z e ) at l e is related to the spreading coefficient S and 
the temperature T throug h 19 ' 21 

4>{l e ) = 5cxT-T w , (7) 

where T w is the wetting temperature, and S is defined by 



S = 7 S 



7si - 7, 



(8) 



where 7 SV and 7 S ] are the substrate-vapor and the 
substrate-liquid surface tensions, the first-order wetting 
transition from the incomplete wetting of a thin liquid 
film with thickness l e (4>(l e ) < = 00)) to the 

complete wetting of a infinite thickness with I = 00 
(</>(Z e ) > = 4> (I = 00)) will occur at the wetting point 
W at T = T w along // = 0" in Fig. [3] Therefore, the 
complete- wetting regime with S > is realized above the 
wetting temperature T > T w . Then, from the Young's 
formula 



7 cos V a = 7 SV - 7si 



we have 



COS^q = 1 



S 



(9) 



(10) 
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FIG. 2: Reduced full potential cj>(l) (Eg. ((23)) ) as a function 
of the reduced film thickness I (Eq. (J24J)) of a completely- 
wettable substrate for various reduced chemical potential ji 
(Eq. (24])). (I) at the prewetting line fi = /z p (/i p = -0.00299), 
(II) below the bulk coexistence and above the prewetting line 
> jti > fip (fi = —0.0015), (III) at the bulk coexistence 
fi = (fi = 0.0000), (IV) above the bulk coexistence fi > 
(fi = +0.0015). At the prewetting chemical potential, the 
thin (with thickness l e ) and thick (with thickness L e ) wet- 
ting film can coexist. The thickness of the thick wetting film 
L e diverges above the bulk coexistence /i > 0. The d = 2 
dimensional droplet height L is determined from the energy 
conservation law (|18p . 
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FIG. 3: Schematic surface phase diagram in the spreading 
coefficient S or the temperature T and chemical potential /i 
plane. The horizontal line /i = corresponds to the bulk co- 
existence line, on which at the wetting point W with 5 = 
(T = T w ) a first-order wetting transition occurs. The prewet- 
ting transition line /i p starts at W and terminates at the 
prewetting critical point PW. Above this prewetting line and 
below the bulk coexistence (0 > fi > /i p ), the stable state is 
a thick film. The thickness diverges at the bulk coexistence 
and this infinitely thick film becomes bulk liquid phase above 
the coexistence (/i > 0). Below the prewetting line (/i < /i p ), 
the stable state is a thin film. A red arrow indicates the route 
in the complete-wetting regime where the isothermal hetero- 
geneous nucleation is expected. 



Therefore the apparent contact angle vanishes (9 a = 0) 
in the complete-wetting regime with S > and, so does 
the free-energy barrier Eq. ([2]). 

Now, we will consider the surface phase diagram in 
the complete- wetting regime (S > 0, T > T w ). In 
this regime, the prewetting transition fi p (T)(< 0) line 
appears below the bulk coexistence line fi = (Fig. [3]) . 
At the prewetting chemical potential \i — /i p , the full 
potential (f)(1) has double-minimum shape with the same 
depth shown as the curve (I) in Fig. [2] which 



cj>{l e ) = 4>{L e ), 



and 



dl 



dct> 



= 0. 



(11) 



(12) 



Then the thin (thickness l e ) and the thick (thickness 
L e ) wetting film can coexist at fi = fi p < 0, which 
plays the role of the bulk coexistence fi = for the 
heterogeneous nucleation in the complete- wetting regime 
(T > T w ). When the chemical potential is increased 
(0 > fi > fip) along the arrow indicated in Fig. [3j the 
thick film with thickness L e becomes stable and the thin 
film with l e becomes metastable above the prewetting 
line (<j)(L e ) < (j>(l e ), curve (II) in Fig. [3]). On the other 
hand, the thin film becomes stable and the thick film 
becomes metastable below the prewetting line (/i < fx p ). 



The metastable thin film loses stability at the upper spin- 
odal h = /j,£ p and the metastable thick film loses stability 
at the lower-spinodal /i = shown in Fig. [3] Finally, 
at and above the bulk coexistence /i > along the arrow 
in Fig. [3J the thickness of the stable thick film diverges 
(curves (III) and (IV) in Fig.©. 



Incidentally, previous authors payed most attention to 
the critical nucleus of wetting transition when [i < 0^1 
and the critical phenomena near the critical point W and 
the prewetting line [ijr^—- We rather pay attention to 
the heterogeneous critical nucleus of bulk phase transi- 
tion in a complete- wetting regime (S > 0, T > T w ) along 
the arrow indicated by "nucleation" in Fig. [3] across the 
coexistence [i > and /i < 0. In an incomplete- wetting 
regime (S < 0, T < T w ) with a finite contact angle 
9 a > 0, a cylindrical and semi-spherical liquid droplet 
(Fig. [IJ can appear on the substrate covered by a thin 
liquid film above the coexistence /i > 0. This is in fact 
the critical nucleus of the heterogeneous nucleation and 
Eq. © of CNT is qualitatively correct. 
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III. CRITICAL NUCLEUS ON A 
COMPLETELY- WETTABLE SUBSTRATE 

A. d = 2-dimensional nucleus 

Based on the phase diagram shown in Fig. [3] and the 
morphology of the potential (f)(1) shown in Fig. [2j we can 
discuss the morphology and the work of formation of the 
critical nucleus on a completely-wettable substrate. Since 
we are most interested in the global character of the het- 
erogeneous nucleation in the complete- wetting regime, we 
will neglect the fluctuation effect that can be important 
near the wetting point W. 

For a d = 2 dimensional cylindrical droplet, the Euler- 
Lagrange equation for the free energy (Eq. ([5])) is simpli- 
fied to2£2S 



dx V(l + ^) 1/2 , 



7^2 



dV 



(1 + '.?) 



3/2 



dl 



(13) 



(14) 



where l xx = d 2 l/dx 2 . Equation (fT4)) could be considered 
as a kind of equation of motion for a classical particle 
moving in a potential —(f>(l). 

Equation (fT4"]) can be integrated once to give 



"7 



2\l/2 



-7 cos 6>(0 = V(l) -fil + C, (15) 



where C is the integration constant and cos9(l) is the 
cosine of the angle 9(1) made between the tangential line 
of the liquid- vapor surface at the height l(x) and the 
substrate^. 

Near the substrate, the liquid-vapor interface of the 
droplet will smoothly connect to the surrounding thin 
liquid film of thickness I = l e with l x = 0, the integration 
constant C in Eq. ([T5]) will be given by C — — 7— V (l e ) + 
fil e , and the liquid vapor interface will be determined 
from 



"7 



(l + ll) 1/2 



(V(l)-V(l e ))-»(l-l e )-j. (16) 



Similarly, at the top of the droplet with a height I = L, we 
have again l x = at / = L (Fig. [2]), and the liquid vapor 
interface will be determined from an equation similar to 
Eq. (jT6]) with l e replaced by L: 



-1 



(V(l)-V(L))-»(l-L)- r (17) 



(l + ll) 1 ' 2 

Since Eqs. (JTH]) and (|17|) must be identical, we have 

V(l e )-fil e = V(L)-(,L, or 4>(l e ) = (j>(L), (18) 

which is similar to the energy conservation law for a clas- 
sical particle whose (pseudo-)cquation of motion is given 



by Eq. (fl"4]) moving in a potential surface —(f)(1). Then, 
the height L of the d = 2 dimensional cylindrical droplet 
can be determined from Eq. (fT5|) . 

Therefore a cylindrical droplet can exist on a 
completely-wettable substrate. The height L of the d — 2 
dimensional droplet is determined from Eq. (fT5]l . Also, 
it exists even in an undersaturated vapor below the bulk 
coexistence and above the prewctting line (0 > fi > fi p ) 
in the thick film phase (Fig. [3J). In this thick film phase, 
the droplet is not the critical nucleus of the bulk phase 
transition, but is in fact the critical nucleus of the thin- 
thick surface phase transition. This critical nucleus is 
expected to transform continuously into the bulk critical 
nucleus above the bulk coexistence (/x+ > > 0) when 
the chemical potential cross the coexistence /j, = 0. 

In order to study the film thickness l e , L e and the 
droplet height L more quantitatively, we use a model 
interface potential 



where lo is a typical thickness of the thin-film, and Vb > 
plays the role of the so-called Hamaker constan t 18 ' 31 ^4 s i v 
of the substrate-liquid- vapor system through 




Volt 



12tt : 



(20) 



and the parameter b plays the role of the tempera- 
ture that controls the transition from incomplete- to 
complete- wetting. Since the two minima of Eq. (|19l) lo- 
cate at I/Iq = 1 and I/Iq = oo, and the spreading coeffi- 
cient S is related to the parameter b through 



S = V(l e ) = Vt 



2-6 



(21) 



the complete- wetting with S > is realized when b < 2. 
We will use the potential parameter b = 1.7, which has 
already been used in Fig. [2J 

From Eq. ([T^)l . the full potential (Eq. ©) can be writ- 
ten as 

0(O = W(i) (22) 
using non-dimensional reduced potential <f> defined by 

(23) 



1 1 + 6 6 

(f)(1) = — = h - ul, 

2/ 2 3Z 3 4l 4 



and 



(24) 



Figure [5] shows the reduced potential <f> (tj in the 

complete- wetting regime for various reduced chemical po- 
tentials fl below the bulk coexistence fi = and above the 
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prewetting /2 p . We set b = 1.7 for which the prcwetting 
chemical potential is given by fx p = —0.00299326. 

Figure 0] shows the reduced stable thick- film thickness 
L e = L e /lo, the metastable thin-film thickness l e = l e /lo, 
and the droplet height L = L/Iq of the d = 2 dimen- 
sional cylindrical nucleus determined from Eq. (fT5|) . The 
height L must be equal to the thick-film thickness L e 
and remains finite at the prewetting line (p = /i p ). We 
observe that those quantities change continuously at the 
bulk coexistence (// = 0) when the chemical potential fi 
is increased from negative (fx < 0, undersaturation) to 
positive (fx > 0, oversaturation) . 




-0.004-0.003-0.002-0.001 0.001 0.002 0.003 0.004 0.005 
Chemical potential p. 

FIG. 4: The stable thick-film thickness L e — L/Iq and the 
metastable thin- film thickness l e = l e /lo and the droplet 
height L = L/Iq of the d — 2 dimensional (2D) cylindrical 
nucleus determined from Eq. (|18[) . The droplet height of the 
d — 3 dimensional (3D) semi-spherical nucleus determined 
from the boundary value problem of the Euler-Lagrange equa- 
tion (|47[) is also shown. 



Using the expansion V (I) — V (L) ~ (I 
L) {dV/dl)\ l=L , we have from Eq. JT7J 



"7 



(1 + ^) 1/2 



dV 
dl 



l=L 



fx )(l - L) - 7 (25) 



near the top of the droplet I ~ L. If the effective chemical 
potential defined by 



(26) 



at the droplet top with height L is positive, then Eq. (|T7|) 
becomes 



"7 



-/Keff (I - L) - 7, 



(1 + ^) 1/2 

whose solution is a semi-circular shape (Fig. [1]) 



off 



x 2 - (i? cff - L) 



(27) 



(28) 



where i? c ff — L is the shift of the base of circular interface. 
The effective radius R c ff is given by the Kelvin-Laplace 
formula^ 



-Reff = 



7 

Meff 



(29) 



In Fig. 03 we show the effective chemical potential 
Aoff = A'cff^o/Uo (Eq. (|26p ) in the complete- wetting 
regime as a function of the chemical potential fx calcu- 
lated from the model potential (|19|) with b = 1.7. The 
effective chemical potential vanishes not at the bulk coex- 
istence /i = 0, but at the prewetting line fi p . The effective 
circular liquid- vapor surface near the top of the nucleus 
can be maintained even in the undcrsaturated vapor with 
Li < as far as (j, e g > (Eq. (|29|l). The effective positive 
Laplace pressure {fi e & > 0) is produced by the repulsive 
surface potential {dV / dl\ l=L < 0) near the top of the 
droplet even though the vapor itself is undersaturated. 
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FIG. 5: The effective chemical potentials fj, e s = Li e glo/Vo for 
the critical droplet in the complete- wetting regime (b = 1.7). 
The effective chemical potential p, e g vanishes not at the bulk 
coexistence /i = 0, but at the prewetting line /i = /2 P below 
the bulk coexistence. 



At the prewetting fx = fx p , the effective chemical po- 
tential [x c ff vanishes 



Moff 



dV 
~dT 



(30) 



l=L 



from Eq. (|12p since L = L e . Hence, the prcwetting line 
[X = /i p acts as an effective or a shifted bulk coexistence 
with fx e g = for the critical nucleus. In fact, Blossey^i 
has noted that the d-dimensional prcwetting line corre- 
sponds to the (d — l)-dimensional bulk coexistence line. 

From Eqs. (|26[) and (f3TJ)) , the effective chemical poten- 
tial [x ff is approximately written as 



Meff 




{fx-fx p ) (31) 
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near the prcwetting line fi — fx p , and the radius R e g 
[Eq. (|29p ] of the curvature at the top of the droplets di- 
verges as 



R 



off 



7 7 
oc 



/'p 



(32) 



at the prewetting line /i = /i p . This divergence has 
already been predicted by Bausch and Blosse y 20 ' 21 us- 
ing a different definition of the droplet radius. We note 
in Fig. [S] that the effective chemical potential /i e g and, 
hence, the effective radius i? e ff (Eq. (|29| ) at the top of 
the droplet change continuously at the bulk coexistence 
(fi = 0). 

The lateral size and the liquid-vapor interface of the 
droplets can be studied only by solving the Eulcr- 
Lagrangc equation (|14[) . which can be done using the 
standard numerical method such as the Rungc-Kutta 
method. To this end, we have to fix the parameter V0/7 
which can be expressed by using the Hamakar constant 
A siv in Eq. $ZU\) and A M by 



1 



where an empirical formula*! 



7 



.4 



lvl 



247TZ) 2 



(33) 



(34) 



with Do = 0.165nm is used. Suppose we tentatively set 
Iq = 2Dq, and using the combining relatio n 1 ^ 31 



Aiv = ~V A n [ V A s S - V A n. 
Ai v i = An, 



we have 




(35) 



(36) 



which will be V0/7 ~ 0.1 — 10 using typical values of A ss 
and Air&. By using the scaled quantities p, e s = /i c ffWK) 
and i? e ff = R e fi/h, Eq. ([29)1 can be written as 



-Rcff 



1 



(Vb/7)Moff' 



(37) 



Therefore, the lateral size of nucleus which is roughly 
determined from R e g is in inverse proportion to V0/7. 

Figure |6] compares numerically determined d = 2 di- 
mensional cylindrical droplet shapes with the ideal semi- 
circular shapes (Eq. ([28])) with the height L and the ef- 
fective radius i? e ff calculated from Eqs. (|18|) and (|29|) 
when V0/7 = 0.5 and 2.0. The droplet shape deviates 
significantly from an ideal circular shape, in particular, 
below the bulk coexistence /j, < 0. The droplet becomes 
flat and its shape becomes pancake-like^ as the chem- 
ical potential is decreased down to the prewetting line 
/i p . The size of the critical pancake is finite even at the 
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(a)The droplet shape when V0/7 = 0.5. 
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(b)The droplet shape when V0/7 = 2.0. 



FIG. 6: The droplet shape numerically determined from 
the Euler-Lagrange equation Eq. (|14|l using the Runge-Kutta 
method (solid curves) compared with the ideal semi-circular 
shape given by Eq. Q280 (broken curves) (I') near the prewet- 
ting line (p = —0.0029), (II) below the bulk coexistence and 
above the prewetting line > n > /i p (p — —0.0015), (III) 
at the bulk coexistence = (p = 0.0000), (IV) above the 
bulk coexistence fj, > (p = +0.0015). Only a right half of 
the droplet is shown. Note the scale of the vertical and the 
horizontal axes. The droplet shape deviates significantly from 
circular shape and becomes "pancake" as the prewetting line 
is approached — > /i p . 



prewetting line (Fig. [5]) even though the effective radius 
i? off diverges (Eq. (j32])). 

It is possible to define an effective contact angle 6 c s 
(sec Fig. |6|) of a droplet on a completely-wettable sub- 
strate by extrapolating the semi-circular shape (Eq. (|2"8)l ) 
down to the surface of the thin film with thickness l e . 
From the geometrical consideration shown in Fig. [TJ we 
find 



cos6> c ff = 1 



L 

Res 



(38) 



7 



from Eq. (|29|) . Therefore, the effective contact angle van- 
ishes (6 B g — > 0) at the prewetting line because L — > L e 
and i? c ff — > oo as /i — > [i p from Eq. (|32|) . 

The work of formation W in Eq. ([T]) can be calculated 
by inserting the droplet profile in Fig. [6] obtained from 
the Euler-Lagrange equation (fT3|) into Eq. ([5]) with d = 2 
and subtracting the contribution from the wetting him 
with thickness L: 



W 



K(i+e) 1/a -i; 

+ (yQ)-fO)-(y(i e )-fii e )]dx, 

which can be rewritten using Eq. (|16p as^i 

7^ 



(39) 



ir 



l=L 



«=«. (1 + ^) 1/2 

2 7 A0 (Z)-A0(J) 2 ) 1/a <fl, 



wherc 



A0(l) = ^(O-^(Ie) = 0(O-^(i). 



(40) 



(41) 



Equation (f40j) is also known as the line tension 16 ^ 11 i 28 . 
Here we interpret this energy as the work of formation of 
the critical nucleus on a completely-wettable substrate. 
Intuitively, the critical nucleus in the complete-wetting 
regime is approximated by a thin fiat disk (cf. Fig. [5]) 
and its free energy is given only by the line tension of its 
perimeter. 

The reduced work of formation W = W/7 from 
Eq. (|40|) as a function of the reduced chemical poten- 
tial jl is shown in Fig. [7(a). In this complete- wetting 
regime, the work of formation W does not vanish and 
does not agree with the prediction of CNT. Furthermore, 
it changes continuously even at the bulk coexistence at 
fi = though the character of the droplet changes from 
the critical nucleus of the prewetting surface phase tran- 
sition below bulk coexistence ^ p < \i < to the critical 
nucleus of the heterogeneous bulk phase transition above 
the bulk coexistence n > 0. Therefore, the nucleation 
rate J given by Eq. ([1]) is expected to change continu- 
ously at the bulk coexistence as weli^. A more detailed 
discussion on the continuity of the work of formation W 
as well as that of the derivatives d n W/dfi n at /1 = will 
be given in the Appendix. 

This work of formation W approaches a finite value 
at the prewetting line fi = /j, p (Fig. [7(a)) due to d = 2 
dimension as the integral (Eq. (|40p ) remains finite. For 
such a flat nucleus (Fig. [6j) with l x <C 1, we can approxi- 
mate (1 + ll) 1/2 -l~l 2 x /2m Eq. (|3H]), and Eq. (gDJ) can 
be approximated hy^ 



W 



l=L 



(2 7 A0 (0) 1/2 dl. 



(42) 



i=i c 



By using the scaled potential in Eq. 

Won hV ) 1/2 l - 



we observe 



(43) 
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(a)The work of formation in the complete-wetting regime. 
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(b)The work of formation in the incomplete-wetting regime 



FIG. 7: (a) The reduced work of formation W = W/'y calcu- 
lated from Eq. (|40[) and the approximate formulas Wcnt = 
Wcnt/7 calculated from Eqs. (|45[) (a) in the complete- 
wetting and (b) in the incomplete-wetting regime. 



Therefore, the reduced work of formation W = W/7 is 
proportional to the parameter \fVopi in Eq. (|33)l (cf. 
two curves in Fig.UJa)). 

Eq. (|40| can also be used to calculate the work of for- 
mation in the incomplete-wetting regime. The result is 
shown in Fig. [71(b) when b = 2.5 and V0/7 = 0.5 and 
2.0. We also show the results using the CNT formula ([2]) 
and ([3]), which will be derived from Eq. ([4"0"]) . Using the 
definition of the local contact angle 9 (I) in Eq. (|16p . we 
can easily transform Eq. (|39|) into simple and compact 
form 



W 



l=L 



smO (I) dl. 



(44) 



i=i c 



By neglecting the surface potential V(l) and using the 
change of variables dl = j sin Odd/ fi from Eq. (fTS"]) , we 
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have 

2 ,fl=*-fl„ 

W C nt = — / sin 2 9d9, 

V Je=e a 

= R 2 cS n{9 a - cos 9 a sin 6 a ) , (45) 

for the incomplete-wetting regime where we have used 
Eq. (|29| with (j, e g = fj,. Equation (|44)) is exactly the work 
of formation in Eq. ^ and ((31) as Whomo = 7ri? 2 ff /^. The 
apparent contact angle # a is calculated from the formula 

cos e a = i + t^—^l (46) 

7 

derived from Eqs. (|15[) by neglecting the surface potential 

The work of formation W calculated from Eq. (|40p is 
always smaller than the work of formation Wont calcu- 
lated from the classical formula (J45j , which is due to the 
presence of the attractive surface potential V(l) which 
enhances the condensation of vapor and lower the bar- 
rier of nuclcation. Both W and Wcnt diverges at the 
bulk coexistence [i = 0. Therefore, CNT is qualitatively 
correct in the incomplete- wetting regime. 

The work of formation in the complete- wetting regime 
in Fig.[7ta) is smaller than that in the incomplete- wetting 
regime in Fig. Efb). However, the former quantity is 
still appreciable and as far as W/k^T > 76 the free en- 
ergy barrier of nucleation should be observable 2 - even in 
the complete- wetting regime. Also, this free energy bar- 
rier and, therefore, the nucleation rate calculated from 
Eq. should change continuously^ at the bulk coexis- 
tence. 



B. d = 3 dimensional nucleus 

The profile of the droplet on the substrate is deter- 
mined by minimizing the free energy (Eq. (|5[)) with re- 
spect to the film thickness I which leads to the Euler- 
Lagrange equation. For a d = 3 dimensional hemispher- 
ical droplet, it is given by2£ 

( d !\ ( lx \ dV 

H^ + *H(i7^J = ^~ M ' (47) 

where l x = dl/dx and x is the coordinate measured from 
the center of the base of the droplet (Fig. [1]). This equa- 
tion cannot be integrated to give Eq. (|T3"|) . as for the 
d = 2 dimensional cylindrical droplet. As a result, it is 
not possible to calculate the height L of the nucleus from 
Eq. (fTB")) . Instead, the height L is determined from the 
solution of Eq. ([47]) that satisfies the boundary condition 
l x = at x = and I = l e at x = oo. 

In Fig. [8] we show numerically determined droplet 
shapes for various [i when V0/7 = 2.0 . In Fig.2]we have 
shown the droplet height L directly determined from the 
numerically determined droplet shape in Fig. [5] as the 



function of the chemical potential \i. By repeating the 
argument from Eq. ([23]) to (j2"9")l , it is possible to approx- 
imate the droplet shape determined from Eq. (|4"T)) near 
the top of the droplet by the semi-spherical shape given 
by Eq. (|28j) . However, the effective radius is now given 
by 

27 

Rcb = — (48) 

Meff 

instead of Eq. (|29|) . Then, it is possible to define the ef- 
fective contact angle 9 c g for a d = 3 dimensional droplets 
on a completely- wettable substrate shown in Fig.[5J This 
effective contact angle is also given by Eq. (|3"5j) and van- 
ishes (9 c g — > 0) at the prewetting line because i? c ff — > 00 
as — > /i p . 



7 























■ (V) \\ \ \(rv)\ (ni)\ 


\<JD 











50 100 150 200 250 300 

Radial distance x 



FIG. 8: The droplet shape numerically determined from the 
Euler-Lagrange equation (|47|l using the Runge-Kutta method 
(solid curves) when V0/7 = 2.0 compared with the ideal 
semi-spherical shape given by Eq. (|28[) (broken curves) for 
(I') near the prewetting line (jl = 0.0027), (II) below the 
bulk coexistence and above the prewetting line > /i > /i p 
{fi = -0.0025), (III) jl = -0.0015, (IV) at the bulk coexis- 
tence fj, = (jl = 0.0000), (V) above the bulk coexistence 
fi > {fi, = +0.0015). The droplet shape cannot be deter- 
mined near the spinodal as its lateral size diverges. Only a 
right half of the droplet is shown. Note the scale of the vertical 
and the horizontal axes. 

In contrast to the d = 2 dimensional case, the lateral 
size of the d = 3 dimensional droplet is larger and di- 
verges at the prewetting line /i = /i p . Using an analogy 
to the classical mechanics, this divergence can be intu- 
itively understandable. Since the Euler-Lagrange equa- 
tion (|47|) can be interpreted as an equation of motion of 
a classical particle^ in a potential field —4>{l) and the 
term (l x /x) plays the role of friction, it takes infinitely 
long "time" x for a classical particle starting from the one 
maximum of the potential —<fi(l) at L = L e to reach the 
another maximum with the same height at l e (sec Fig. [3] 
at the prewetting when j2 p = —0.00299). Therefore the 
"lateral size" x of the d = 3 dimensional critical droplet 
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diverges at the prewetting line. For a d = 2 dimensional 
droplet, since there is no friction term and the energy 
conservation is assured from Eq. (|T8|) . the lateral size of 
the droplet remains finite at the prewetting. 

Since we can not integrate Eq. (|T7|) analytically for a 
d = 3 dimensional droplet, we cannot use Eq. (|40[) to 
calculate the work of formation W. Instead, we have 
to resort to a direct numerical integration using Eq. (|5|). 
Figure H] shows the reduced work of formation W = W/ 7 
of a d = 3 dimensional critical nucleus as a function of 
the reduced chemical potential fl. 




-0.004-0.003-0.002-0.001 0.001 0.002 0.003 0.004 0.005 
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FIG. 9: The reduced work of formation W = W/~f calculated 
from Eq. ((5]) by using the droplet profile l(x) numerically de- 
termined from Eq. (|47|) when V0/7 = 2.0 as a function of 
the reduced chemical potential jx. In contrast to the d — 2 
dimensional nucleus, the work of formation W of the d = 3 
dimensional nucleus diverges at the prewetting fi p . 

Similar to the d = 2 dimensional critical nucleus, the 
work of formation IF is a continuous function of the 
chemical potential \x at the bulk coexistence. However, 
in contrast to the d — 2 dimensional nucleus in Fig. E£a), 
the work of formation W of the d — 3 dimensional nu- 
cleus diverges at the prewetting /i p in Fig. [5J 

IV. CONCLUSION 

In this paper, we have used the interface displace- 
ment model to study the heterogeneous nucleation on a 
completely-wettable planar substrate where the first or- 
der incomplete to complete wetting surface phase tran- 
sition takes place. The prediction of the classical nu- 
cleation theory (CNT) is critically tested. It is found 
that the classical picture breaks down on a completely- 
wettable substrate, where CNT predicts that the critical 
nucleus and the nucleation barrier are expected to van- 
ish since the apparent contact angle vanishes. In fact, 
both the critical nucleus of the heterogeneous bulk phase 
transition in the oversaturated vapor and the critical nu- 
cleus of the surface thin-thick prewetting transition in the 



undcrsaturation vapor can exist and transform continu- 
ously at the bulk coexistence on a completely-wettable 
substrate. 

Therefore, the critical nucleus as well as the free-energy 
barrier of nucleation exist on a completely-wettable sub- 
strate. Furthermore, the nucleus exists even under the 
undcrsaturatcd vapor. The undcrsaturated vapor turns 
to an effectively oversaturated vapor due to the effec- 
tive interface potential of the substrate which prefers in- 
finitely thick wetting layer. The Laplace pressure (or the 
effective chemical potential fi c g in Eq. (f2U)) ) that is neces- 
sary to produce semi-circular liquid- vapor interface of the 
droplet becomes positive by the interface potential even 
under the undcrsaturatcd vapor (Fig. [3]). Then, the bulk 
coexistence does not play any role to the critical droplet 
in the complete- wetting regime. Instead, the prewetting 
line plays the role of the coexistence in the complete- 
wetting regime. Various properties of the critical nucleus 
changes continuously as functions of the chemical poten- 
tial. Therefore the nucleation rate (Eq. (JTJ)) is expected 
to changes continuously at the bulk coexistence. Our 
result support the conclusion reached from a numerical 
simulations in Ising system by Sear—. He simply stated 
that this small nucleus docs not "know" whether it will 
grow to form a wetting layer of finite thickness or a bulk 
phase of divergent thickness. Our re-examination of var- 
ious properties of the critical nucleus in the complete- 
wetting regime supports his conclusion and, hopefully, 
will rekindle the interest on heterogeneous nucleation on 
a completely-wettable substrate. 
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Appendix A: Continuity of the work of formation W 
and its derivatives d n W/djj, n a.t the coexistence n = 

Since, the full potential <p(l, /1) is non-singular at fi = 
as shown in Fig. 2J all quantities including the solution 
I of the Euler-Lagrange equations (fTB"]) and (|4"T)) is ex- 
pected to be non-singular at \x = 0. Therefore, the work 
of formation W calculated from Eqs. and (|40[) is a 
continuous function of the chemical potential (i at the 
coexistence fi = as shown in Figs.JTJa) and [5] 

For a d = 2 dimensional cylindrical droplet, we can 
check not only the continuity of W but also that of its 
derivative d n W/dfj, n directly from Eq. (|3TJ|). First, the 
work of formation W in Eq. ([40|) is a continuous function 
of the chemical potential \x at fi = as indicated by 
Fig. EJa) since the thin-film thickness l e (fi) determined 
from dcj)/dl — 0, the droplet height L (fi) determined from 
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Eq. Q8|). and the integrand of Eq. (HO 



from Eq. (|A2I) . we find 



to(i,/i)= (2 7 A0(Z, A/ )-A0(Z,/.) : 



1/2 



(Al) 



as a function of the chemical potential \x are all non- 
singular at /I = (see Eq. ([6]) and Fig. SJ. 

Also, the derivatives dW/d/i can be continuous (dif- 
fcrentiablc) at /i = as long as the derivatives dl e /dfi, 
dL/dfi, and dw(l,(i) jd\i exist and are all non-singular 
at fi = 0. For example, from d(j> (l e ) /dl = 0, we find 



dV (le 
dl 



/'■ 



By differentiating Eq. (| A2|) . we obtain 



die 
(/// 



d 2 V(l ) 
dl 2 



(A2) 



(A3) 



where Iq = l e (fi = 0). Therefore dl e j d\i exist at /i = 
and is non-singular as far as the right-hand side 
of Eq. (|A3|) is finite. By successively differentiating 
Eq. (|A2[) . we can easily prove that the n-the derivative of 
thin-film thickness d n l e /d[j, n exists and is given generally 
by 



d n L 



dfi n 



In 



d 2 V(l ) d n+1 V(l Q ) 



di 2 



(A4) 



for n > 1, where /„ is a rational function. Then, n-the 
derivative d n l e /dfi n will be non-singular at /j, = as far 
as all derivatives d k V (Iq) /dl k (k = 1, 2, . . . , n + 1) are 
finite. 

A similar argument can be applied to the droplet 
height L (fi) determined from Eq. (fl~8f . Again by differ- 
entiating Eq. (HU) and using Eq. ([A"3j) and (Iq) /dl = 



dfx 



(Lq-Iq) 



dl 



(A5) 



exists and is non-singular at (i = 0, where Lo = 
L (/i = 0). Again, we can easily prove by mathemati- 
cal induction that the n-th derivative of droplet height 
d n L/dfx n is given generally by 



d n L 



= 9n ( I 



d 2 V(l Q ) d n V(l ) 



dl 2 ' ■ ' ■ ' dl n 
dV(Lo) d n V(L ) 



dl 



dl n 



(A6) 



for n > 2, where g n is a rational function. Therefore, the 
derivatives d n L/dfj, n exists and is non-singular at u = 
if all the arguments of g n in Eq. (|A6[) arc non-singular. 



Since all the partial derivatives dw k (l,n) /d/i k (k = 
1, . . . , n) arc non-singular at /i = 0, the derivative 
d n W/dfi n will be continuous at /i = as far as all the 
derivatives d k l e /d[i k and d k L/d/j, k (k = l,...,n), and, 
therefore, those of the effective potential d k V(lo)/dl k 
(k = 2,...,n+ 1) at l and d k V (L ) /dl k (k = l,...,n) 
at Lq are non-singular. In particular, since our model 
potential Eq. ([T9| is infinitely diffcrcntiable at Iq and Loj 
all derivatives d n l e /d^ n and d n L/d/j, n are expected to 
be continuous at = (Fig. 0]). Therefore, not only 
the work of formation W but also its all the derivatives 
d n W/dii n will be continuous at /i = 0. 
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